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Abstract: In the present letter we find that Starobinsky's inflationary solution is also 
valid in the Dvali-Gabadadze-Porrati (DGP) model where a 3-brane is embedded in 5- 
dimensional Minkowski bulk. We show that such a solution is typically not supported 
by the Self-Accelerated branch of the model, giving therefore a natural selection of the 
conventional branch of solutions. In the absence of brane induced Einstein-Hilbert term 
the SA branch is always selected out. We then study the linearized modes around all such 
de Sitter brane solutions finding perturbative stability for a range of parameters of the 
brane QFT. 
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1. Introduction 

The DGP model of brane induced gravity a model for modification of gravity at large 
distances (cosmology at late times) provides an appealing alternative explanation for the 
current accelerated expansion of the Universe @, i- 

In this article we study the very early Universe in brane induced gravity. Inspired 
by the seminal work of Starobinsky Q we consider the localized matter as being the high 
energy limit of an asymptotically free QFT (such as QCD) and search for inflationary 
de Sitter solutions (see for a more recent discussion about Starobinsky's solution). A 
thorough study of quantum effects of brane conformal field theories in the physics of the 
early (brane) Universe, has already been done for several other models of brane world: it 
includes the RS model as well as other variants that involve the inclusion of higher 
curvature bulk terms, a richer bulk matter content, as well as super symmetric setups and 
generalizations to dS bulks and FRW brane cosmology 

The main motivation of the present study of inflation in brane-induced gravity mod- 
els, relies on the observation that if the DGP model successfully incorporated inflation 
then, interestingly enough, it would mimic the whole cosmological paradigm: an initial 
meta-stable inflationary solution that would decay into "regular" FRW cosmology. In the 
high energy limit, when the brane curvature satisfies R > M^^r^, the brane stress ten- 
sor is well approximated by the trace anomaly-generated stress tensor quadratic in brane 
curvatures ^ acting as a source for the junction condition associated to the 4d modified 
Einstein equations of motion. Such a junction condition will, as usual, lead to the brane 
Friedman equation. 

We analyze the cosmology resulting from the model, with and without the inclusion 
of an induced non-conformal EH term on the brane. We then proceed to study the per- 
turbative spectrum and find the conditions on the parameters of the brane quantum field 
theory such that there is a period of inflation that is long enough. 
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In the original DGP model there is a branch of "self-accelerated" solutions [|| in which 
late time acceleration could be driven solely by the effect of the modified gravitation without 
the need for any energy density. This branch of solutions, whose consistency with data has 
been extensively studied (see, e.g., |9|), appears to be unstable both perturbatively [10, |ll| 
(see, however, the possible caveats to this conclusion |T^) and non-perturbatively (with 
hints in the exact solutions 14 1 and a recent discussion in |]l5[ - see also, p6[). A welcome 
feature that results from our approach is that once the quantum effects that drive inflation 
are included the conventional (as opposed to the self-accelerated) branch of solutions of the 
model is naturally selected, if the induced Planck mass is small. On the other hand if the 
induced Planck mass is large the self-accelerated branch results perturbatively unstable. 



2. The Setup 

We consider a five-dimensional action which includes the Einstein-Hilbert term Seh and 
a bulk matter action Sb which we will assume to give rise to a perfect fluid type of stress 
tensor Sf^ = —2g^^^5SB / Sg^^ = <iiag{—pB,PB,PB,PB,Py)- We will be interested in 
studying brane-cosmological bulk solutions of the form 

ds^ = -N\t, y)dT^ + A\t, y)dnl + B\t, y)dy^ (2.1) 

where y is the coordinate orthogonal to the hypersurface {y = 0} that we will refer to 
as "the brane" and = Qrnn{x"^)dx"^dx^ is a unit radius maximally-symmetric space, 
characterizing the brane spatial section. We will assume that some matter is confined on 
the brane. In particular, we will consider the backreaction caused by the renormalized 
stress tensor of a set of massless conformal quantum fields on an isotropic homogeneous 
(brane) universe, along the lines of what originally done by Starobinsky 0]. Such massless 
quantum fields are to be understood as the high energy limit of the (massive) fields of 
an asymptotically free QFT. The result will be twofold: on the one hand we generalize 
Starobinsky's idea to the DGP brane world model. On the other hand we study brane- 
induced gravity (DGP model Q) in case when the matter stress tensor on the brane is 
dominated by the trace-anomaly-induced one. 

It was shown [jl^ that, assuming there is no bulk-brane energy exchange, the rr and 
yy bulk equations of motion take the simple form ^ 

F' = -^A'A'pB (2.2) 
F = '^AA% (2.3) 



where 



nt,y) = [^) -1^1 -kA^, (2.4) 




^We denote ' — dr and ' — dy 
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whereas the ry equation of motion reads 



N' A A' 



N I A 
A N 



. 



(2.5) 



Equation ( ^.2| ) can be integrated and yields 

^2 



F 



6 



■PbA'' - C , 



(2.6) 



where C is an integration constant that was found |18| to be related to the bulk Weyl 
tensor. Hence, 



AB 



AN 



+ 



k 



C 



A^ A^ 6 



■PB 



(2.7) 



We will be interested in the case of an empty, Minkowski bulk for which = C = so 
that F = Q and the leftover Einstein's equation is also automatically satisfied and (| 
reduces to 



Hence, 



i- 

\AB 



A' 
AB 



AN 



+ 



k 

l2 



(2.8) 



AN 



+ 



A^ 



(2. 



which projected upon the brane becomes 



a' 
ah 



a\ k 

- +— ' 
a J 



where 



a(r) = ^(r,0) , 
6(r) = i?(r,0) , 
n(r) = N{t, 0) = 1 



(2.10) 



(2.11) 
(2.12) 
(2.13) 



The l.h.s. of ( p.lO| ) is related to the extrinsic curvature on the brane that is defined by 



where 



M 



dX 



M 



(2.14) 

(2.15) 
(2.16) 
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are respectively the tangent vectors characterizing the brane embedding in the bulk (we 
choose a static gauge) and the (unit) normal vector to the brane 



nM = (0,0, B) 



(2.17) 



so that 



jy-m 



n 
J 

-I 
ab 



+ 



±1 , 



(2.18) 
(2.19) 



is the expression for the extrinsic curvature tensor, where we have made use of ( p. 10] ). The 
Israel junction condition links the latter to the stress tensor of the matter localized on the 
brane 

1 



2M3^^^ ' 



(2.20) 



where the Ihs is evaluated at y = 0+ and 7^jy is the induced metric on the brane, namely. 



-dr^ + airfdnl = -i^y{x)dx^'dx'' . 



(2.21) 



Assuming the brane matter to be a perfect fluid, from the spatial components of the 
extrinsic curvature we get 



a\ k 



1 



6M3 



P ■ 



(2.22) 



If the localized matter is conformal the (v.e.v. of the) energy momentum tensor reduces 
to the one generated by the trace anomaly ||^, pOj 



= a 



'6 



-2V^,V^i? + 2gf,^aR + 2RR^^ - ^g^^^R^^ = aSj,^ + /5e^^ , (2.23) 



where a = A;2/6!(27r)^ and /3 = k^/6l{2TT)'^ and all tensors are calculated according to the 
induced metric, and yield 



k 

a J 



1 



2M3 



a 



a\ a a\ a 
- +2- - +2- 

a a \a a 



a\ k 

2 



-a] +^ + U +^ 



(2.24) 



that is the brane Friedman equation which can be written as 



2M3 



VVaa VVa v'l , kVa 



(2.25) 
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by means of the transformation v{a) = [aa)"^ + ka?' . One possible set of solutions of the 
latter equation can be achieved by setting v = H'^a^ (with constant), upon which the 
term in the round parenthesis vanishes (lower "indices" a denote derivatives w.r.t. a); the 
leftover terms uniquely fix the constant, 

{aaf + ka^ = H'^a^ , 

aH^ = 2M^e , (2.26) 

and select the allowed branch, that is e = +1 for k2 > and e = — 1 in the opposite case. 
In other words de Sitter type of solution 



a(r) 



oo eyip{Ht), k = 

H-^cosh{Ht), k = +l (2.27) 
H'^smh{Ht), k = -l 



holds in our brane world setup. However, since ^ 

k2 = Ns + llNp + 62Nv > , (2.28) 

only the conventional branch (e = +1) is allowed for such a solution. In fact, a negative 
contribution to ^2 can be obtained from quantum fields that have a wrong sign in the 
kinetic term, in other words ghosts: since the trace anomaly in four dimensions comes from 
a triangle diagram describing the correlation function of three stress tensors, changing the 
sign of the action for the quantum fields changes that of the stress tensor and in turn that 
of the anomaly. The previous argument clearly raises an issue concerning the (quantum) 
stability of the branch e = —1, in the present setup. Such a branch corresponds to the 
so-called self-accelerated branch of the original DGP model Q. ^ Hence, the argument 
we gave above would a priori seem to confirm the more pessimistic points of view. Note 
however that such an argument does not imply that the self-accelerated branch is altogether 
excluded by ( p. 24 ); in fact the second term on r.h.s. of (2.24) is not positive definite - it is 



related to the trivial anomaly - and thus it might allow for solutions even for e = — 1. 

Before concluding this section let us also point out that in the above solution the 
inflation rate H is set in by the fundamental scale M and depends upon 4(i physics only 
through the numerical combination k2- In the (DGP) brane world scenario the fundamental 
scale is taken to be much smaller that the 4d Planck mass: therefore the curvature scale 
R ~ is also much smaller than Mp and thus lower than M^yrp, so that the validity of 
the approximation ( p. 23 ) appears not to be guaranteed. However, we will see in the next 



^Ns is the number of real scalar fields, Nf the number of Dirac fermions and Nv the number of vector 
fields (21). 

^Let us stress that here we refer to e = +1 (e — —1) as the conventional (self-accelerated) branch just 
to make contact with analogous definitions given in DGP cosmology Q where the classification is related 



to the sign of the total energy density as in (2.22): our reference to such names is just nomenclature and 
is not linked to the properties of the solutions. In fact, in our case both dS branches are, strictly speaking, 
self-accelerated meaning that the acceleration is geometrical and not driven by sources. For the sake of 
clarity in fact, the "true" conventional branch would be flat 4d Minkowski that is obviously a solution 



of (2.22) 
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section that the inclusion of an induced non-conformal Einstein-Hilbert term on the brane 
enhances the value of the de Sitter scale. A brane-induced E-H term is also a necessary 
ingredient in order to recover 4d gravity between brane massive sources. 



2.1 Adding a non-conformal Einstein-Hilbert term on the brane 

In this subsection we consider the inclusion of a non-conformal induced EH term on the 
brane. Such an induced term is part of the effective action of matter fields coupled to 
gravity [0. However it was shown that a brane EH term can be produced also at the 
classical level on a tensionful brane if the bulk action includes higher curvature terms |^2[ . 

The presence of an induced E-H term on the brane does not necessarily mean that 
the inflationary phase is no longer there. In fact as we now show we can have a steeper 
inflationary phase whose scale is set in by a combination of both M and induced Planck 
mass Mp. Note in fact that the anomaly generated stress tensor, "ameliorated" with the 
localized Einstein tensor coming from the induced EH term, is a good approximation for 
"large" brane curvature, Mp > R > Mq^j,. The first inequivalence guarantees that brane 
loop corrections are under control. Note also that, although bulk loops are controlled by 
the fundamental scale M, the expansion is under control there as well since the bulk is (a 
patch of) flat Minkowski space. The aforementioned inflationary phase could be unstable 
as pointed out by Starobinsky and Vilenkin ^] and the universe would thus finally decay 
to "regular" FRW cosmology. We discuss about stability in the next sections. 

In this section we thus consider the inclusion of an Einstein-Hilbert term on the brane 



Ml 



/ 



d Xy/—gRi 



so that ( p. 261) becomes 



erric + H 



a 



Ml 



(2.29) 



(2.30) 



with nic = 2M^/Ml being the crossover scale of the DGP model. Denoting z = H/M\, 



with A = (2a) equation ( 2.30| ) can be rewritten as a cubic equation 



2A 



MpV 
M J 



z 



4e = 



(2.31) 



that is of the form that allows for direct Cardano's solution 

1/3 



MA 



2e + 



+ 



2e 



2A (Mp^ 



1/3 



. (2.32) 



For e = +1 there is one real positive solution regardless of the value of the quantity 
X{Mp/Mf. 

For e = —1 and < \{Mp /M)^ < 3 • 2~^/^, there is no positive solution whereas 
for \{Mp/M)^ > 3 • 2"^/'^ there are two positive solutions. Therefore, if the induced 
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Planck mass is large enough, both branches would seem to allow expanding de Sitter solu- 
tions ( ^.271) . In fact, in the limit where Mp ^> M one recovers Starobinsky's expansion rate 
Hs ~ Mpl\fdL. Note also that the expansion rate for the conventional (self-accelerated) 
branch corresponds to > 1 {^^r ^ other words a phenomenological selection 

"criterion" that, a priori, seems to pick the conventional branch is the highest dS expan- 
sion rate. However, in such a case higher order corrections (in H/ A) to the quantum stress 
tensor are not guaranteed to be under control, since the induced Plack mass is related to 
the UV cutoff of the localized theory as Mp ~ A^A^, with N being (roughly) the total 
number of fields in the theory [N < k2) and thus Hs ^ A. An hypothetical way out to 
such a problem might be the inclusion of bulk higher curvature terms that could allow to 
lower the dS scale. Let us also remind the reader that, in order for quantum gravitational 
corrections to be under control, Hs < Mp, the total number of fields must be extremely 
large such that A:2 > 6! ^ 



3. Perturbations of the de Sitter brane solution 

In this section we consider small perturbations around the de Sitter brane solution intro- 
duced previously. Let us first point out that the stress tensor 0^j,, as defined in ( 2.23| ), is 



not conserved in a generic background and should be supplemented by a term -R^^C^ojy/j, 
with being the Weyl tensor. However, linear perturbations of such a term around 

a maximally-symmetric background, vanish: we can thus safely omit the latter in what 
follows. 

The bulk equations and junction conditions are given by: 

Rmn = , (3.1) 
G,. - me {K,, - g,.K) = ^ + ^ , (3.2) 

where we added to the r.h.s. of (|3.2| ) an extra subleading T^y component. 

We want to solve for the perturbations around the /c = de Sitter background solution 
(see the appendix) in a convenient gauge, namely, 

g^J.u = N'^ivhfiu + h^u , 9yy = '^ 9^ = ^^ Vh^y = Vi,h. (3.3) 

where N{y) = 1 — eH\y\, has curvatures R^j.^ = 3H'^^f_iu and R = 12H^ and e = ±1 
(e = -|-1 is the conventional branch: \y\ < H~^). 
Let us first gather some results: 

5R^.u = \y^Nvh - V + V " H'^li^u'h , (3.4) 
5R = -3H^h , (3.5) 



* A severe parametrical constraint comes from the observational limit on the amplitude of long- wavelength 
gravitational waves, requiring fcs 3> ^2 that appears quite unnatural. However notice that fcs is the 
coefficient of the trivial part of the anomaly and is thus arbitrary. 
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Thus, 

5G^. = (nV - V^V,/i) - 2/7^ V + ^H'^f,,h , (3.6) 



From the yy equation we get ^^'^5Kfj_y = —eHh/2 fl^ (cfr. eqs. (15-16)), which combined 



with the background equation ( 2.30 ) gives 

5 {-mc{K^y - gfj^uK)) = -\mcdyh^y - ^ ^ H'^-f^yh + 4(1 - a)H'^'h^y , (3.7) 



where a = aH^/Mp and, for future reference, /? = PH'^/Mp. Therefore, the junction 
conditions ( ^I^ ) and its trace yield: 

-l i^cdy + (1 - 2a - 4/3)n - AH\l - |a - 2/?)^ V = 



^ (r^. - ^7/..^) - ^ ((1 - 2a - 6/5)V^V, - (1 - 2a - 2(5)H^-i^,) h , (3. 



Mp 
1 



-ei/m, - ai/^ + /3nj /i = . (3.9) 

Note that for non vanishing (3 the trace h becomes a propagating degree of freedom. 

A way towards the solution of this system is through the shift to transverse traceless 
metric perturbations h^y by introducing the brane bending mode (p{x) [^] 

2N 

V = V + ' (3-10) 

with = v^V/,±7^^//2^ 



The system (|3.8| ), (^) and the bulk iiv equations becomes: 

-\ (rricdy + (1 - 2a - 4/?)n - AH^l - - 2/3) ) V 



J/ 

where we use 



^ (t^. - ^7/..r) + {-em, - 2H{a + 4/3)) n+ - 2|:n+ Qy. , (3.11) 

HT 

{eHm, + 2ai72 - 2/3n) Q^^ = ^ , (3.12) 

{N^d^ - O) V = ' (3-13) 



and the identity 



O = 4if2 - □ , (3.14) 
Q = -4H^ - □ , (3.15) 



On+ = U-,Q<I) . (3.16) 
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valid for any scalar <p. Combining ( 3.11 ) and ( |3.13| ) into a single equation gives: 



rricdy - 
-45(y) 



+ O - 26{y) {{2a + 4/3)n - 2H\3a + 4/3)) 



1 



1 



Tf^u - ^If^uT j - {em, + 2H{a + 4/?)) n+ 99 - 2-^n+ 



(3.17) 



which by making use of the condition ( 2. 301) derived from the background equations, i.e., 
erric = {a — 1)H, can be solved to yield 

h 

2 

where, 



S{0) (t^, - ^Tm-t) - {em, + 2H{a + 4/?)) n+ - 2-|n+ Q(^) (3.18) 



(1 -a){z- H^q{z)) - (q + 4/3) {z - 2H^) ' 



5(^) 



Q(z) = -(l + eJl + 4-^ 



(3.19) 
(3.20) 



We recover /i^j^ via ( |3.10D with ip given by the solution of (|3.12| ) 



Notice that we have correlated the two possible solutions (3.20) with the two branches 
of the model via the e parameter. This implies a choice of boundary conditions at y = 
00: had these been uncorrelated, non-normalizable modes would be contributing to the 
amplitudes implying the presence of sources far away into the bulk. 

Having the solution for the perturbations ( p. 18] ) we now concentrate on the one gravi- 
ton exchange amplitude between brane localized sources T and T' (we follow the approach 



presented in IQ; see [24] for a pedagogical introduction). To this end, it is useful to decom- 
pose the corresponding conserved energy momentum tensors Tf^i, and T'^'^ into transverse 
traceless, longitudinal and trace parts: 



T, 



flU 



T, 



TT 



1 



1 



1 



(3.21) 



where Py,y = V^V,y — ^"ff^uO (for which the identity P^uf{Q,)(t) = f{0)P^y4) holds for any 
smooth / and scalar (f)). 

We obtain the following expressions: 

1 „,„.„TT 



1 



1 



S^HMj, \Q 2(3Q- 



T, 



(3.22) 
(3.23) 



where ^ = 1 — 3a — 8/3. 

Next, let us introduce the Lichnerowicz operator A with the following properties 



(A - 4:H^)T^'^ 



{A-4H^)P^,<P = P^,Q<P 



(3.24) 
(3.25) 
(3.26) 
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This gives, 



5(0)TjJ = S{A - 4H^)T^, - \s{Qh^,T - -V^V 



1, 
3 



T'i^-S{0)T^J^ = T'^^SiA - 4ij2)r^^ + Lr Q - 3^ S{Q)T + V^(- • •) . 
We compute the ampUtude, 



A 



(3.27) 



(3.28) 



obtaining: 



A 



Ml 



2f3N 



jj2 \ jY 



(3.29) 



Let us focus on the conventional branch (e = +1). In this case, there could be poles 
at Q = 2iJ^ (vanishing denominator of S{Q)), at Q = and at 2/3Q = ^H'^. We will also 
discuss below the interpretation of the other zero of the denominator of S{Q). 

The Q = 2ij2 (□ = -6i?^ = -2A) pole, together with the one at A = GiJ^ corresponds 
to a massless spin-2 state on the de Sitter background (with the correct tensorial structure). 

The would be Q = (□ = —4H^) pole has vanishing residue (at that point H^S{0) + 
N/^ = 0). 

To go further into the analysis, let us find another expression for the amplitude. We 
define a complex variable w = —z — H^/A such that S{w) = S{—w — has a branch 

cut in the positive real w axis. We get the alternative form 



ziri Jy s — w 



-Sis) 



lim (s — wq)S{s) 



(3.30) 



Wq — W s^wo 

via a contour integral T which goes around the cut just below and above the positive w 
axis and closes at infinity through a counterclockwise circle. The residue of the pole in the 
first Riemann sheet, wq, plus the jump in the imaginary part across the cut give 

'•^ ds 



S{w) = 



Wq — W 

3 



+ 



/ 



W 



2 - 5a - 12/3 ' 

^ (1 - a^H^s + ((1 -2a- 4/3)s + 3i?2(i _ ^ _ g^^^s 



(3.31) 
(3.32) 
(3.33) 



or equivalently, 

SiQ = -D- 4ij2) = 



-□-2A 



+ 



dm- 



2m 



□ -2A 



p{rr? -^H'^/A) . (3.34) 
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The simple pole at w = wq = —9H^/4: in ( |3.31 ) corresponds to the massless state (□ + 2A 



in ( ^.34 ), see equivalence with expressions in [^]) while the integral has contributions 
from the continuum of massive KK modes < w < oo, i.e., 9-ff^/4 < □ + 2A < oo with 



positive definite spectral density ( 3.33| ) as long as a > 1 (a resonance representing the other 



zero in the denominator of S" - see ( |3.19|) ). With these conventions, the massless state has 
positive norm as long as c is positive, which for the regime of interest for us, namely, a ^ 1, 
implies /3 ^ —1/4. 

The 2/3Q = £^H^ pole has positive residue (corresponding to an extra scalar state with 
positive norm) as long as ^ is positive, which for a ^ 1 coincides with the positivity of 
norm of the massless state. The squared mass (eigenvalue of □ + 2H^) of this extra state 
is 

ml = -2H^ - . (3.35) 



By recalling the relation to the parameters k2 ( 2.28| ) and ^3 in Q 

we see that a modest negative fcs partially compensating ^2 in ( |3.36| ) leads to a theory 
free of ghost-like instabilities around a de Sitter background. Furthermore, in the original 
Starobinsky setup a future instability is found due to a light tachyon. This slow 

instability which could be responsible for the exit from a long inflationary phase persists 
in our setup. Keeping ^ > and focusing in the regime a ^ 1, /3 ^ —1/4 gives a tachyonic 
mass ^ —2H^ to the extra scalar state. 

On the self-accelerated branch (e = —1) the pole w = wq sits in the first Riemann 
sheet provided one chooses the opposite prescription for yj—w |12|. With this choice we 



obtain that the spectral density has the same form as in (g) though now is negative 



definite, since a < 1: it can be interpreted as due to a tower of negative-norm KK states. 



4. Conclusions 



We considered the high energy limit of asymptotically-free brane matter in Minkowski 
bulk and studied Starobinsky's trace-anomaly driven infiationary solutions. We studied 
perturbations about such de Sitter solutions and found that the system could display 
no instabilities. In particular, in the conventional branch there is a range of parameters 
(dictated by the number of scalar, vector and fermionic components) of the brane QFT 
such that there are no ghost-like excitations but a slow instability to end the inflationary 
phase. 

In absence of any induced brane EH term, we found that only the conventional branch 
allows the aforementioned solution, becoming naturally selected. In the presence of an 
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induced EH term the phenomenon is more complex: for small induced Planck mass, Mp < 
Md^/^, only the conventional branch is allowed, whereas for big induced Planck mass 
Mp > Md^/^ both branches are allowed. In particular in the limit Mp ^> M both 
branches converge to the Starobinsky's solution. For the self-accelerated branch we have 
H < Hs that corresponds to a < 1 (a = 1 being Starobinsky's solution), whereas for the 
conventional branch we have H ^ Hg and a ^ 1. However, the self-accelerated branch 
turns out to be unstable regardless of the value of Mp. 

Acknowledgments 

It is a pleasure to thank F. Bastianelli, G. Gabadadze, N. Kaloper, S. Odintsov and L. Sorbo 
for very useful discussions and the referee for very helpful criticisms. We are grateful to the 
organizers of the Simons Workshop 2006 held at SUNY at Stony Brook, where this work 
was initiated, for their hospitality. AI was supported by DOE Grant DE-FG03-91ER40674. 
OC is grateful to the Department of Physics of UC Davis and to the Institute of Cosmology 
and Gravitation of the University of Portsmouth for support and hospitality while parts of 
this work were completed. 



A. Bulk solutions 

We briefly sketch the derivation of the bulk solutions corresponding to the dS cosmologies 
studied above. As shown by ||l^, ^ the assumption of zero bulk-brane energy exchange 
(cfr. eq. (^)) allows to write 

4 = . (A.l) 

Thus setting the brane "clock" such that A^(t, 0) = 1, and using the bulk equation of 
motion (^) we get 



A{T,y) = a(r) - elylVk+l^ (A.2) 

N{T,y) = l-e\y\^=^ (A.3) 

and therefore, using the explicit solutions described in the main part of the paper we obtain 

A{T,y) = a{T){l-eH\y\) (A.4) 

N{y) = 1 - eH\y\ (A.5) 

and thus 

ds^ = N{yfdf + dy^ = N{yf ( -dr^ + a^ir) dnU + dy^ , (A.6) 



with \y\ < H ^ for the conventional branch e = +1 and \y\ < oo for the self- accelerated 
branch e = — 1. In conformal coordinated the latter reads 

ds^ = e-'"\'\ {-dr"^ + (r) d^l + dz^^ (A.7) 

with z e M. 



- 12 - 



References 

[1] G. R. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. B 485, 208 (2000) 

[arXiv:hep-th/0005016]; G. R. Dvali and G. Gabadadze, Phys. Rev. D 63, 065007 (2001) 
[arXiv:hep-th/0008054]. 

[2] C. Deffayet, Phys. Lett. B 502, 199 (2001) [arXiv:hep-th/0010186]. 

[3] C. DcfTayct, G. R. Dvali and G. Gabadadze, Phys. Rev. D 65, 044023 (2002) 

[astro-ph/0105068]. 

[4] A. A. Starobinsky, Phys. Lett. B 91, 99 (1980). 

[5] J. C. Fabris, A. M. Pelinson and I. L. Shapiro, Nucl. Phys. B 597, 539 (2001) [Erratum-ibid. 
B 602, 644 (2001)] [arXiv:hep-th/0009197]. 

[6] S. W. Hawking, T. Hertog and H. S. Reall, Phys. Rev. D 62, 043501 (2000) 

[arXiv:hep-th/0003052]; S. Nojiri and S. D. Odintsov, Phys. Lett. B 484, 119 (2000) 
[arXiv:hep-th/0004097]. 

[7] S. Nojiri, S. D. Odintsov and S. Zcrbini, Phys. Rev. D 62, 064006 (2000) 

[arXiv:hep-th/0001192]; S. Nojiri, O. Obregon and S. D. Odintsov, Phys. Rev. D 62, 104003 
(2000) [arXiv:hep-th/0005127]; S. Nojiri and S. D. Odintsov, JHEP 0007, 049 (2000) 
[arXiv:hep-th/0006232]; S. Nojiri and S. D. Odintsov, Phys. Rev. D 64, 023502 (2001) 
[arXiv:hcp-th/0102032]; S. Nojiri and S. D. Odintsov, JHEP 0112, 033 (2001) 
[arXiv:hep-th/0107134]; K. Koyama and J. Soda, JHEP 0105, 027 (2001) 
[arXiv:hep-th/0101164]. 

[8] A. Vilenkin, Phys. Rev. D 32, 2511 (1985). 

[9] C. Deffayet, S. J. Landau, J. Raux, M. Zaldarriaga and P. Astier, Phys. Rev. D 66, 024019 
(2002) [arXiv:astro-ph/0201164]; A. Lue, R. Scoccimarro and G. Starkman, Phys. Rev. D 69, 

044005 (2004) [arXiv:astro-ph/0307034]: Phys. Rev. D 69, 124015 (2004) 
[arXiv:astro-ph/0401515]; J. S. Alcaniz, D. Jain and A. Dev, Phys. Rev. D 66, 067301 (2002) 
[arXiv:astro-ph/0206448]; E. V. Linder, Phys. Rev. D 72, 043529 (2005) 
[arXiv:astro-ph/0507263]; L. Knox, Y. S. Song and J. A. Tyson, arXiv:astro-ph/0503644. 
G. Dvah, A. Gruzinov and M. Zaldarriaga, Phys. Rev. D 68, 024012 (2003) 
[arXiv:hep-ph/0212069]; M. Fairbairn and A. Goobar, Phys. Lett. B 642, 432 (2006) 
[arXiv:astro-ph/0511029]; K. Koyama and R. Maartens, JCAP 0601, 016 (2006) 
[arXiv:astro-ph/0511634]; I. Sawicki and S. M. Carroll, arXiv:astro-ph/0510364; R. Maartens 
and E. Majcrotto, Phys. Rev. D 74, 023004 (2006) [arXiv:astro-ph/0603353]. 

[10] M. A. Luty, M. Porrati and R. Rattazzi, JHEP 0309, 029 (2003) [arXiv:hep-th/0303116]; 
A. Nicolis and R. Rattazzi, JHEP 0406, 059 (2004) [arXiv:hep-th/0404159]; K. Koyama, 
Phys. Rev. D 72, 123511 (2005) [arXiv:hep-th/0503191]; D. Gorbunov, K. Koyama and 
S. Sibiryakov, Phys. Rev. D 73, 044016 (2006) [arXiv:hcp-th/0512097]; K. Koyama and 
F. P. Silva, Phys. Rev. D 75, 084040 (2007) [arXiv:hcp-th/0702169]. 

[11] C. Charmousis, R. Gregory, N. Kaloper and A. Padilla, JHEP 0610, 066 (2006) 
[arXiv:hep-th/0604086]; 

[12] C. Deffayet, G. Gabadadze and A. Iglesias, JCAP 0608, 012 (2006) [arXiv:hep-th/0607099]; 

[13] N. Kaloper, Phys. Rev. Lett. 94, 181601 (2005) [Erratum-ibid. 95, 059901 (2005)] 

[arXiv:hep-th/0501028]; Phys. Rev. D 71, 086003 (2005) [Erratum-ibid. D 71, 129905 (2005)] 
[arXiv:hep-th/0502035] . 



-13- 



[14] G. Dvali, G. Gabadadze, O. Pujolas and R. Rahman, Phys. Rev. D 75, 124013 (2007) 
[arXiv:hep-th/0612016]. 

[15] R. Gregory, N. Kaloper, R. C. Myers and A. Padilla, JHEP 0710, 069 (2007) 

[arXiv:0707.2666 [hcp-th]]. 

[16] K. Izumi, K. Koyama, O. Pujolas and T. Tanaka, Phys. Rev. D 76, 104041 (2007) 
[arXiv:0706.1980 [hep-th]]. 

[17] P. Binetruy, C. Deffayet, U. Ellwanger and D. Langlois, Phys. Lett. B 477, 285 (2000) 
[arXiv:hep-th/9910219]. 

[18] T. Shiromizu, K. i. Maeda and M. Sasaki, Phys. Rev. D 62, 024012 (2000) 

[arXiv:gr-qc/9910076]; S. Mukohyama, T. Shiromizu and K. i. Maeda, Phys. Rev. D 62, 
024028 (2000) [Erratum-ibid. D 63, 029901 (2001)] [arXiv:hep-th/9912287]. 

[19] P. C. W. Davics, S. A. Fulhng, S. M. Christensen and T. S. Bunch, Annals Phys. 109, 108 
(1977); P. C. W. Davics, Phys. Lett. B 68, 402 (1977). 

[20] S. W. Hawking, T. Hertog and H. S. Reall, Phys. Rev. D 63, 083504 (2001) 

[arXiv:hep-th/0010232]; S. W. Hawking and T. Hertog, Phys. Rev. D 65, 103515 (2002) 
[arXiv:hep-th/0107088]. 

[21] M. J. Duff, Nucl. Phys. B 125, 334 (1977); F. BastianeUi and P. van Nieuwenhuizen, Nucl. 
Phys. B 389, 53 (1993) [arXiv:hep-th/9208059]. 

[22] O. Corradini, A. Iglesias, Z. Kakushadze and P. Langfelder, Phys. Lett. B 521, 96 (2001) 

[arXiv:hep-th/0108055]. 

[23] J. Garriga and T. Tanaka, Phys. Rev. Lett. 84, 2778 (2000) [arXiv:hep-th/9911055]. 
[24] G. Gabadadze, arXiv:hep-th/0408118. 

[25] M. Porrati, Phys. Lett. B 498, 92 (2001) [arXiv:hep-th/0011152]. 



-14- 



